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Abstract 

The possibility of composite systems arising out of a point charge interacting with a 
Nielsen-Olesen vortex in 2+1-dimensions is investigated. It is shown that classical bounded 
orbits are possible for certain ranges of parameters. Long lived metastable states are 
shown to exist, in a semi-classical approach, from the study of the effective potential. 
Loss of self-adjointness of the Hamiltonian and its subsequent self-adjoint extension in 
some cases leads to bound states. 
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1 Introduction 



2+l-Dimensional physics has come of age in the last couple of decades. Specifically with 
the advent of anyons [p]] , (excitations having arbitrary spin and statistics), and its posited 
presence ranging from diverse condensed matter systems such as high T c superconductiv- 
ity, fractional quantum Hall effect to more exotic scenarios in the high energy regime such 
as processes in the presence of cosmic strings etc., 2+1-dimensions has ceased to be just 
a laboratory for testing ideas conceived for the "physical" 3+1-dimensions. The present 
paper deals with studying the energy spectrum of a quantum (charged) particle in the 
presence of a Nielsen- Olesen magnetic vortex |2| , in 2+1-dimensions. 

The motivation is twofold. The physical existence of magnetic vortex lines in type II 
superconductors []3| (in an external magnetic field) makes the study of a vortex-particle 
system interesting. Also according to the present lore and as envisaged by Wilczek [p, 
a possible realisation of anyons in nature is some sort of a composite object, consisting 
of a magnetic flux tube (of the "fictitious" Chern Simons gauge field) attached to a bose 
or a fermi particle. Albeit the crucial difference in nature of the two magnetic fields, the 
existence of (quantum) bound states, (as we demonstrate here), should lend credibility to 
the hypothetical anyon structure. 

Let us elaborate a little on our system, our modes of analysis as well as the results that 
we have obtained. The system is that of a non-relativistic point charge in the presence of 
a vortex, the latter being the Nielsen- Olesen vortex solution in the Abelian Higgs Model 
(AHM) H ||. We restrict ourselves to minimal gauge invariant coupling between the 
point charge and the gauge field. 

We start with a thorough discussion on the classical aspect of the problem, which 
turns out to be quite tricky. We show that classical bounded trajectories || of the particle 
are allowed for certain ranges of the parameters of the model. This analysis is important 
since it hints at the possible existence of quantum bound states, which brings us to the 
Schrodinger equation problem of the charge in presence of a classical vortex potential. 
One would be too optimistic to conclude that upon quantization, the closed orbits would 
correspond to bound states. Apparantly this is not the case. Essentially this is because 
the Hamiltonian in question being positive definite, negative energy bound states can not 
appear. Still we will show that the effective potential energy profile is such that it allows 
metastable states of quite long life time. Previous works in a simlar vein are those of 
a point particle in the presence of a 't Hooft-Polyakov monopole in 3+1-dimensions. 

However, bound states in the strict sense do exist, for certain values of effective angular 
momentum at least. This is related to the loss of self adjointness of the Hamiltonian, 
due to the singular vortex potential, and its subsequent self adjoint extension, ]7|, ||, which 
allows only the above bound states. An early work in this connection is 0. 

The paper is organised as follows: in section 2, we introduce the vortex potential and 
the detailed classical dynamics. In section 3, we discuss the metastable quantum bound 
states. Section 4 deals with the problem of self adjointness and the ensuing true bound 
states. We conclude the paper in section 5. 

2 Classical dynamics 

The Nielsen-Olesen vortex solutions || have provided a deep connection between a rel- 
ativistic field theory, (the AHM), and the system of type II superconductors. If the 
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superconducting material is invariant along the ^-direction, the system is essentially a 
two-dimensional one, in the x — y-plane In the former case vortices appear in the 
two space dimensional slice of the magnetic flux lines extended along the third direction. 
Hence one can study the vortex solutions in the 2+1-dimensional AHM. 

Let us start by introducing the vortex solutions of the AHM [fj, [|]. The (c = 1) 
Lagrangian is 

Lahm = --F^F, V + i(ZV>rP M 0) - ^(1 I 2 ~F 2 ) 2 , (1) 

where F^ u = d^A u — d u A^, D^<p = d^cjy — ieA^cf). The problem has a cylindrical symmetry 
and the asymptotic solution, (i.e., the vortex), in polar (r, 6) coordinates is 

A r = A = 0, Ag{r) » — + — , r oo, (2) 

er 



<P = e in \F + f3e-^ Fr ). (3) 

We will be confined to the n = 1 or single vortex line or 'flux-tube' sector, since numerical 
work has indicated it to be stable, a and (3 are constants and a is of the order of ^sE, 

In this work we only consider the particle to be interacting with the gauge potential A^. 
Obviously, other interactions such as Yukawa or vector couplings with can be studied 
as well. Henceforth we consider the gauge field to be external. Clearly this particular 
solution of A^ in (0) generates no electric field. The magnetic field is 

B = -2yr5 (2) (r) + a(— - e F)^—=- (4) 

e 2r 



The matter Lagrangian with U(l) gauge invariant coupling is 

Lmatter = \Mr 2 + ^Mf 2 ^ + (5) 

where j 7 * is the conserved particle current. We denote time and space derivatives of O by 
O and O' respectively. Writing explicitly the potential, we obtain 

1 1 -1 ae~ eFr 

L m atter = -Mr 2 + -MrV + gr9(— + — ), (6) 

2 2 er 



Note that the term qr0(±) = s 6> is a total time derivative and does not influence the 

& \ er ' e 

classical equations of motion. It is dropped in the present classical analysis. The Lorentz 
equation for the particle is 

MI' = ge ij X j B, (7) 

where M and g are mass and charge of the particle and B is without the 5-function term. 
In polar coordinates we have 

X = re r + r6e e , (8) 

X = (r-r6 2 )e r + (r9 + 2f6)e e , (9) 
and the Lorentz equation in component form is 

M(r - r6 2 ) = gr9B, M(r9 + 2r9) = -grB. (10) 
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Note that 2 has the dimension of h. In the quantised version the particle charge g becomes 
related to e by g = k(eh) where k is some number. The 9 equation following from ([]) is 
given by 

Mr 2 6 + ag^e~ eFr = L, (11) 

where L is a constant. Thus a generalised angular momentum is conserved. Borrowing 
from the classical terminology, this is one of the integrals of motion, the other being the 
conserved energy, since a static background interaction is chosen. The r-equation provides 
the energy, when we substitute 9 from (|TT|) 



d 1 ,L — ag\fre t '' ' 
dr^2M 



Mr = --^[— ( "~^ V ' C ) 2 ], (12) 



and the conserved energy E is 



E = -Mr- + J_ { L-a9^e-^ = 1 J_ 

2 2M V r ; 2 2M V ; V ; 

This is just the sum of kinetic and effective potential energy. 

To proceed further, we must ascertain first that the potenial energy function V(r) in 



13|) does have a profile which is able to sustain bounded classical motion of the particle. 
Clearly V(r) is positive infinity at r = and goes to zero at large r. We want to argue 
that between r = and r = oo, there is one minimum and one maximum, which will give 
the well. If we express the potential energy V(r) as V(r) = [T(r)] 2 , then 

V{r)' = 2T(r)T(r)\ V{rf = 2(T') 2 + ITT" . (14) 

Clearly V = has two roots, T = and X" = 0. It is easy to see that the root T = 
constitutes the minimum. To show that under some restrictions, the other root X" = 
does represent a maximum, V" = 2TT" has to be negative. We evaluate V" explicitly 
and use X' = to express it in the following form 



Hence the allowed range of r is such that either | > eFr > or eFr > |. 
Thus a potential well is formed with the minimum ar r\ obeying 

L = ag^F ie - eFr \ (16) 

and the maximum at r2 obeying 

L = ag^r- 2 e- eFr2 (^ + eFr 2 ), (17) 

provided r 2 > r\ and T\ is within the allowed range. Obviously the minimum value of 
Vir) is zero, the expression being positive definite. 

Let us now turn to numerical results. Since the equations for the minimum and 
maximum, (fL6|), fllTD , involve transcendental functions, we have given graphical solutions 
in Fig. 1 with the potential well strucure in Fig. 2, for the particular solution chosen in 
Fig. 1. From Fig. 1 one can see that there is a lower bound for the gradiant of the straight 
line (II), below which there is no solution for t\. However, for this extreme value of the 
parameter, solution for r 2 is such that r 2 < rj.. On the other hand, there are solutions for 
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ri and r 2 , with r 2 > r 1; which is required for well formation, such that r\ is at least less 
than -i, that is the relevant dimensionless variable ei*Vi < 1. In Fig. 1, (II) and (III) 
give the values eFr\ = 0.08, eFr 2 = 0.21. In Fig. 2, we have plotted the dimensionless 
variables V(r)/ vs ' e -^ r - The well is formed within the range of the vortex, where 
the exponential damping can compete with the ^j^s centrifugal term. But in this region 
the asymptotic solution may not be too reliable. 

We now give a rough idea of how the system might look close to the origin. One can 
have solutions of and <fi field equations, (obtained from (1)), such that 



A = A r = 0, A e = aeF 2 r + (3e 3 F 4 r 3 , = CeF 2 r + pe 3 F 4 r 3 , 



with a = (3 = 24(1 A + 4A) , C 2 = -3^x7, V = "if, A = ^. The field equations are 
satisfied to 0((eFr) n ), n > 4. Comparing with [[|, we find that near the origin, <fr goes 
to zero correctly. The magnetic field is 

B = — + 0(r 2 ). (18) 

e 



B will behave properly near origin if the 0(r 2 ) term opposes the constant value, which is 
true for 4r > t- 

Now, the L matter and the energy are 

Lmatter = \mt 2 + ^Mr 2 2 + gr6(aeF 2 r + 0(r 3 )), (19) 

E = \mt 2 + ^(L - ageF 2 r 2 ) 2 . (20) 

Apart from a constant shift of La ^ F this potential energy function is just that of a two 

dimensional harmonic oscillator, with the minimum value of zero at r± = \J L/ (ageF 2 ). 
Hence we note that an approximate solution of the gauge field, valid at short distance, 
provides in the leading order, a two dimensional circular oscillator potential, which indeed 
has stable classical orbits. These two limiting results, that is gauge fields at large and 
small r, do indicate that the system in question is classically stable. 



3 Quantum dynamics: metastable states 

Let us notice at the outset that the Hamiltonian considered by us is of the form 

H=^[p- 9 A(r)} 2 , (21) 

where in our case, A(r) is the vortex potential 0. However, without going into details, 
for a purely magnetic field, if is a non-negative, hermitian operator. We consider the 
Hilbert space to be composed of wave functions which are square integrable on the plane 
and regular at the origin. This means that the energy spectrum is positive and purely 
contineous, contrary to the nature of bound states. We will discuss later in more detail 
the Schrodinger equation analysis which also points to the same conclusion as above. For 
the present, arguing from a physical point of view, a potential well which is everywhere 
positive and goes to zero asymptotically can not produce a bound state, (for special cases 
contrary to this conclusion, see Simon in 0), simply because the particle can tunnel out 
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and live outside the well. This example clearly shows that closed classical trajectories, 
which are present here, are in no way a sufficiency condition for quantum bound states. 
After the above compelling arguments against the existence of bound states in the strict 
sense, (that is states having negative energy), let us look at the potential profile in Fig. 
2 more carefully. We immedietly notice that although the well is quite sharp, (since the 
terms responsible are exponential), the subsequent fall off towards zero is very gradual, 
(thanks to the inverse power law). Obviously this will reduce the tunnelling amplitude 
to a large extent and so the states having positive energies that are well inside the "well" 
can not live outside the well due to the greater potential energy. Thus the well quite 
efficiently traps the particle states and these are termed by us as metastable states or 
resonances. Even though the square integrable eigen functions are disallowed, in the spirit 
of semiclassical quantization, we generate harmonic oscillator states about the potential 
minimum. This will give the spectrum of metastable states. 

It is now straightforward to perform the oscillator approximation near r\. While 
considering the asymptotic solutions, we write r = r% + x in V(r) of (|i~3|) and keep upto 
0(x 2 ) terms and get 

V(n +x)= 2M ^ + X? \L - ag^r 1+ x)e-^^, (22) 
- -^[ag^e-^i^- - eF) 2 ]x 2 . (23) 



2Mr\ v v 2r 
Identifying this energy as \kx 2 we find 

/T j_ _ eF|= jBM (24) 

V M M y/n 2ri M K J 

Note that this frequency corresponds precisely to the Larmor frequency for a charged 
particle in a magnetic field. 

Hence the metastable energy spectrum consists of the levels 

E n = {n + -)hw. (25) 

(The anharmonic corrections can be included as well). 

Near the minimum, the energy of the system in oscillator approximation becomes 

E = -Mx 2 + -kx 2 = —p 2 + -kx 2 . (26) 
2 2 2AT 2 v ; 



Alternatively the Bohr-Sommerfield condition can be used to obtain the spectrum of (|2q), 
which is an equation of ellipse in phase space. Thus / p x dx is the area of the ellipse, which 
is a multiple of h. One can also compute the phase integral J^p r dr where a and b are 
the classical turning points obtained by solving (p0| ) when Mr = 0. However near the 
potential minimum, 

t p r dr = f V2MEdr = V2ME, (b - a) w x. (27) 
All these will give rise to a spectrum similar as that of fl25|). 
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In case of the solutions valid near the origin, once again expanding near r\ and writing 
V{r\ + x) — \x 2 , we find 

™ = \l^= g — (28) 
|M M V ; 

with B computed in flI8|). This is also of the same form as (23). Here we can give the 
exact energy levels or wave functions for the point particle. 

Indeed, one can do a detailed Bohr-Sommerfield quantisation, by introducing radial 
and angular quantum numbers for the r and 9 phase space integrals. The standard way is 
to replace 4 by -4 (using the 6 equation of (|TT[) ) in the r equation obtained from (|), in 
order to obtain an equation for the orbit . The energy E appears as a parameter. Using 
the orbit equation and the definition of the quantum numbers, the energy is expressed in 
terms of the quantum numbers. However, due to the complicated nature of the equations 
of motion, we have adopted a more pedestrian approach to obtain the metastable state 
energy spectrum in a harmonic oscillator approximation about the potential minimum. 



4 Quantum dynamics:true bound states 

Let us start this section with some preliminaries, which are distinct for 2+1-dimensions 
from the conventional 3+1-dimensions, only in some technical details. For a generic central 
potential V(r), the complete set of eigenfunctions are of the form ip(r, 8) = f(r)exp(U9), 
with I an integer for singlevaluedness, where f(r) satisfied the radial equation 

d 2 Id I 2 2M 
dr z r dr r z n 



with /(0) = |10|. In our particular case, with the explicit form of V(r) from (|13|), we 



have to solve the radial equation 

[f* + l T~ «" + ^ eFr f ~ k 2 )}h(r) = 0, (30) 
dr z r dr r eyr 

where L = I - f , leZ and k 2 = 2ME > with f L (0) = 0. Note that here the total time 
dervative term in (|6]) has not been dropped, so that the full expression of B in @) is being 
considered. 

One should not confuse the addition of - in the angular momentum (which came from 
the -9 term mentioned above in the Lagrangian), with the ageold recipe of changing 
particle statistics in the quantum version. In external field problem angular momentum 
is quantised in integers even though the expression differs from the canonical one [|TT|]. 
The effect of total derivative terms in quantum mechanics is briefly discussed in 

In our case, with positive eF, for large r, (that is away from the origin), the exponential 
terms die out quickly leaving only the inverse power law falloff, which controls the large r 
behaviour. Now, if the Hilbert space is that of regular square integrable wave functions on 
the plane, then the spectrum is purely continuous with k > 0. In fact the eigenfunctions 
asymptotically behave as 

Mkr)~CMkr), (31) 

where C\ is a constant |TIJ. The above analysis rigorously shows that bound states in 
general are not possible. 
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However, one must be more careful in imposing boundary conditions |L2) on the wave 
functions due to the presence of the 5-function term in the magnetic field. This points 
towards a contact interaction which visciates the self adjointness of the Hamiltonian at 
the origin. The remedy is to look for self adjoint extensions (if present) of the Hamiltonian 
and this can lead to bound states |J for some restricted set of parameters. 

Generally we tend to overlook the difference between Hermitian and self adjoint op- 
erators in quantum physics, but there is a difference in the structure of their respective 
domains, which crucially governs the dynamics in many cases of physical interest. For 
unbounded operators, (that do not have bounded expectation values, such as energy, mo- 
mentum, angular momentum and position), one has to specify the domain of the operator 
as well as its acton on the domain. Without geting involved in too many technical details, 
an operator T is called Hermitian if T C T*, that is D(T) C D(T*) and Tip = T*0 for 
all <peD(T). Here the adjoint T* is defined as (Tip, 0) = (ip,T*(p) for all ip,<peD(T) with 
D(T) denoting the domain of T. T is self-adjoint if T = T* that is iff T is Hermitian 
and D(T) = D(T*). Only the self-adjoint operators can be exponentiated to give one 
parameter unitary groups which governs the dynamics in quantum mechanics. 

The basic criterion of self-adjointness is the following: Suppose T is a self-adjoint 
operator and a exists such that T*0 = i<p. Then also Tip = i<p and 

-z(0, <P) = (i<f>, <P) = (T0, 0) = (0, T*0) = (0, T0) = i{<j>, 0) 

which implies = and this is true for T*<p = —i<p as well. Hence for self-adjointness, it 
must be ensured that ker(T* ±i) = or D(T*) = D(T). However, if T is not self-adjoint 
on a domain, but has n± independent solutions T*0j = ±i0j for some <pieD(T*) then, only 
if n + = n_ = n, one is allowed to make an extension of T to r, where D(T) C D(t) and 
rip = Tip for all 4>eD(T)), such that r is self-adjoint. Basically one modifies the domain of 
T to r such that D(r) = D(t*) and this also ensures that ker(r ± i) = 0. The vectors 0j 
generate the deficiency subspace and n± are referred as the deficiency indices. Note that 
for simple differential operators, as in the present case, T* is the same as T, but acting 
on D(T*), which in general is different from D(T). 

Our task is now to express the Hamiltonian in a suitable form T and find independent, 
normalisable solutions of 

T0^ = ±z0i, (32) 

Since we are looking at regions close to the origin, the exponential term is dropped, 
compared to more singular centrifugal term. Replacing f(r) and kr by and p in fl2T3|), 
we obtain 

d 2 u p p(p + 1) 

— 2" + (1 a — > u p = °> ( 33 ) 



dp 2 p 2 



where p(p+ 1) = — (t — L ). In general the solutions of the Schrodinger equation as well 
as its first order partial derivatives will be continuous, uniform and bounded functions 
over all space, including p = 0. Here we also have /(0) = 0. This induces the following 
natural boundary conditions on u p as well; 

uJO) = = 0. (34) 

Let us see somewhat heuristically why the singularity problem at r = affects only some 
of the p states. The general solution of fl33|) near the origin is 



u p w ApP +1 + Bp~ p . (35) 
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For p different from or -1, to maintain regularity at p = 0, either A or B has to vanish 
and u p (0) = 0. Thus the singularity problem does not matter for these states. But for p 
equal to or -1, Uq or tt_i can be finite but non zero at the origin. For them the imposition 
of the boundary condition u(0) = clashes with self adjoint property of the Hamiltonian. 
This problem and the subsequent solution is elaborated below. 



Let us factorise ( J33D 



d p+l d p + 1 

)u p = v p , (— H )v p = -Up, (36) 



dp p dp p 

or 



— + -)u p = v p , (- )v p = -Up, (37) 



dp ,dp 

^ + -)Up = Vp, {- 

dp p dp p 

For p taking values or 1, the pairs of equations become identical and reduce to 






(38) 

The deficiency subspace is generated by a pair of normalised vectors (i.e. n + = n_ = 1) 

T X ± = ±m, X± = e' r ( ±- J ■ ( 39 ) 
Now the extension of T is r where T = r with domain 

D(t) = Uj + P X+ + IPX- I ^eD(T),(3eC}, (40) 
and the operation of r on D(r) is 

r(V> + Px+ + IPX-) = T^ + iPx+ ~ vyfc-, (41) 

where | 7 | is an isometry that maps x+ X-- We parametrise 7 = e ia and there are 
diferent extentions for different a. In this extended domain the deficiency indices vanish, 

n±(r a ) = 0. (42) 

This is the self-adjoint extension explained before. The vital role played by the contact 
term in the magnetic field in (|4]), which gave rise to the s term in L becomes manifest 
only now. Note that in the expression p(p + 1) = — (i — (I — f ) 2 ), since I is an integer, 
the LHS can vanish, (i.e., p = 0, —1), only if the - term cancells | in the RHS. 

From a slightly different angle, let us now see how the boundary conditions clash 
with the self-adjointness of the Hamiltonian at the origin. Indeed, this is probably an 
easier way to understand the connection between deficiency subspaces and the bound 
state spectrum in question. From the definition of the self adjoint operator, it follows 
that the Hamiltonian is self adjoint if 

(0,TV>)-(T0,VO = O 

with (f) = ( Ul j eD(T*) and <fi = ^ U<2 ^ eD{T). Using the natural boundary conditions 

1*2(0) = 1*2(0) = and W2(°o) = ^2(00) = 0, (since these specify D{T)), the above 
condition reduces to 

ui*(0)« 2 (0)-« 1 *(0)v 2 (0))=0 (43) 
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Thus the above equation vanishes, independent of any bounday condition on <j>, which is 
in the domain ofT*. Thus D(T*) > D(T) since it consists of similar vectors as in D(T), 
but without the boundary condition. This makes T not self adjoint in this domain. The 
cure is to moderate the too strong boundary condition on D(T) by eg., 

u 2 (0) + akv 2 {0) = u 2 (0) + em' 2 (0) = 0, (44) 

where a is an arbitrary real parameter. This is the extension r. This ensures identical 
boundary conditions on ip and <j>, and hence the domains D(t*) and D(t) are identical, a 



and a of (|4l|), fl44]) , are related in the following way, 

la , 
a = -cot-. (45) 

Now we can get the bound state spectrum easily. For either p = or p = -1 we have 

u = kv, v' = —ku, (46) 
We also have from the extended domain 

u(0) + au'(0) = 0. (47) 

Using the above equations we get 

- rv'(0) + akv(0) = 0. (48) 
k 

Let us consider some form of u and v which vanish at r = oo. 

u = Pe~ 7r , v = Qe~ ir . (49) 
P, Q and 7 are constants. Solving the above set of equations we find 

*» = -! 

a 2 

, and so the energy is 

E = -ih- (50 > 

Note that we used the two component factorised form from (|33D since the (matrix) 
differential operator in fl5S|) is Hermitian in its domain, although individually i-4- is Her- 
mitian. However, since ^ is Hermitian, we can use (|33|) directly as well. Choosing p = 
or -1 with 



u = Ae^ r , u(0) +au'(0) = 

we again obtain k 2 = —-%. The difference in these two formulations is nontrivially mani- 
fested in the deficiency index analysis where the domains of these operators are different. 

Thus (BD) is the cherished expression for the bound state spectrum. These states 
correspond to the parameter values such that L = I — 2 = | with integer values for I. So 
far we only know that a is real. There are an infinite possibility of different extensions 
for distinct values of a and still further conditions are required to uniquely specify a. For 
a zero energy bound state, we require a = 00 and u'(0) = but there is no restriction on 
u(0). 



10 



5 Conclusion 



We have considered a system of point charge, interacting with a Nielsen-Olesen vortex 
of the abelian Higgs model, in 2+1-dimensions. We have shown that classical bounded 
orbits are allowed, for certain restricted range of parameters. 

Next we quantise the particle motion, treating the potential as external. From gen- 
eral arguments it is established that in general negative energy bound states are not 
present. However, the nature of the effective potential energy allows metastable states of 
considerable lifetime. 

On the other hand, presence of the contact interaction, (due to the vortex), makes the 
analysis subtle for certain angular momenta. We have shown that at least for some values 
of the parameters, at the origin the self adjoint property of the Hamiltonian is lost, which 
is restored by the self adjoint extensions, and in this new domain, bound states appear. 
There is a one parameter arbitraryness in the bound state energy spectrum which can be 
removed by introducing further physical input. 

Recently we have come across the papers fl3| where the self-adjoint extensions of the 
Dirac Hamiltonian have been studied, in the context of Aharanov-Bohm effect and in the 
presence of a magnetic vortex. 
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Figure Captions 



Fig. 1 



I — > e x vs. x, II — > ax vs. x, III — > + 1) 2 t>s. x, where x = 2eFr, a = |(^) 2 = 

7.4 



Fig. 2 



V(x) = - v/^oie- 11 ) 2 with a = 7.4 
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